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Abstract 

First-principles calculations were performed to investigate the effect of Pd concentration and oxygen vacancies on 
the stability of Pd at LaFe0 3 surfaces. We found a much stronger tendency of Pd to segregate by taking the 
aggregation of Pd at LaFe^Pd^Os-y surfaces into consideration, resulting in a pair of Pd-Pd around a vacancy. 
Moreover, we predicted that one oxygen-vacancy-containing Fe0 2 -terminated surfaces would be stable at high 
temperatures by comparing the stability of LaFe^Pd/Vy surfaces, which further supports our previous conclusion 
that a Pd-containing perovskite catalyst should be calcined at 1,073 K or higher temperatures in air to enhance the 
segregation of Pd in the vicinity of surfaces to rapidly transform the Pd catalyst from oxidized to reduced states on 
the perovskite support. 
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Background 

A three-way catalyst simultaneously transforms toxic 
exhaust emissions from motor vehicles into harmless 
gases. However, the sintering problem, i.e., the growth 
and agglomeration of precious metal particles on con- 
ventional catalysts during vehicle use dramatically de- 
grades catalytic activity, and large amounts of precious 
metals are required to retain the activity of catalysts after 
long periods of use. Thus, intelligent catalysts have 
attracted worldwide attention due to their greatly im- 
proved durability as a result of the self-regenerative 
function of precious metal nanoparticles [1-3]. It has 
been confirmed that the activity of catalysts can be 
preserved, and the amount of precious metals that are 
required can be reduced by 70% to 90% [4,5]. The 
self-regenerative function, which can be explained 
as resulting from the transformation of the state of 
precious metals (Pd, Pt, and Rh) that reversibly move 
into and out of the LaFe^JVi^Os perovskite lattice, 
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significantly suppresses the growth of precious metals 
during the use of catalysts. 

Thus far, many experiments have been devoted to 
research on the state of Pd in perovskite in redox pro- 
cesses. Uenishi et al. [6] investigated the superior start- 
up activity of LaFePdO^ at low temperatures (from 100°C 
to 400°C) using X-ray spectroscopic techniques under 
the practical conditions where they controlled automo- 
tive emissions. They found the Pd° phase partially segre- 
gated outside the surface even at low temperatures; thus, 
the segregation of Pd° under a reductive atmosphere 
induced the start-up activity of LaFePdO^. Eyssler et al. 
found a high concentration of Pd distributed on the 
LaFe0 3 (LFO) surface that contributed to high methane 
combustion due to the formation of PdO in which Pd 2+ 
was in square planar coordination. Additionally, two Pd 
species (Pd 2+ at the surface and Pd 3+ in a solid solution) 
were found to be generated in further calcination. Pd 2+ 
and Pd 3+ could be transformed in equilibrium under 
thermal treatment conditions [7,8]. More recently, Eyssler 
et al. studied the state of Pd in different B-site substitutions 
and compared the effect of catalytic activities on methane 
combustion. A well-dispersed octahedral Pd-O species was 
found for Fe- and Co B- site cations, and PdO particles 
were on the LaMn0 3 surface [9]. Above all, related investi- 
gations have become more important as the activity of 
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catalysts strongly depends on the state of the precipi- 
tated Pd. 

Hamada et al. [10] more recently found in a density 
functional theory (DFT) investigation that oxygen vacan- 
cies (V Q s) that formed near the LFO surface promoted 
the segregation of Pd. They also proposed a scenario 
that perovskite containing precious metal is calcined 
during the catalyst preparation step at 1,073 K for 2 h in 
air, and then V Q s are produced that enhance Pd segrega- 
tion, resulting in a LaPdC^ layer that eventually forms 
close to the surface. The LaPdC^ layer in the vicinity 
of the surface promotes efficient switching between Pd 
metal particles under reductive conditions and the 
dissolved state of Pd in the LaFe!_^Pd^0 3 perovskite lat- 
tice under oxidative conditions. Therefore, the LaPdC^ 
layers formed in the vicinity of the oxide surface play a 
key role in the self-regenerative function. Almost simul- 
taneously, transmission electron microscopy observa- 
tions [11] of atomic-scale processes in Pd-LFO catalysts 
have demonstrated that redox reactions between the 
formation of Pd particles on the Pd-LFO surface under 
reducing conditions and the dissolution of Pd particles 
into LFO under oxidizing conditions take place in 
spatially-limited areas, especially in the proximity of 
oxide surfaces, indicating strong interactions between Pd 
and oxide surfaces. Katz's results also provided strong 
support for the mechanism proposed by Hamada et al. 
However, the stability of the LaPdOs^ layer and the 
mechanism for Pd leaving the LaPd03_ 7 layer have not 
been discussed in detail. The interaction between Pd 
atoms in the perovskite host is especially important 
considering the possibility of nanoscale spinodal decom- 
position as pointed out by Kizaki et al. [12]. Therefore, 
we systemically studied the relative stability of the 
Pd m V 0w -containing surfaces (m =1 and 2 and n =0, 1, 
and 2) in our present work to investigate possible 
phases appearing in steps to prepare catalysts at high 
temperature in air. 

Methods 

Model and computation 

We have calculated the lattice constants [13] of LFO 
and the segregation tendency of Pd at two terminations 
of the perovskite surfaces with and without V Q by using 
state [14,15] and quantum ESPRESSO (QE) [16] codes. 
We found that both state and QE codes yielded the 
similar bulk lattice constants and caused the segregation 
behavior of Pd, which was a strong indication that both 
codes could admirably describe the properties of Pd 
incorporated in the LaFe^Pd^Os^ surfaces. Here, we 
employed the state code to do the first-principles 
calculations. The ion-electron interactions were de- 
scribed using ultrasoft pseudopotentials [17], and the 
exchange and correlation potential was represented by 



a generalized gradient approximation (GGA) in the 
Perdew-Burke-Ernzerhof formula [18]. DFT calculations 
with Hubbard correction (DFT+U) are known to correct 
the bandgap and magnetic moment in local- density ap- 
proximation and generalized gradient approximation cal- 
culations. This method can yield reasonable agreement 
with the experimental results. We omitted DFT+U from 
this work because Hamada et al. verified [10] that elec- 
tronic structures with DFT+U are qualitatively the same 
as those in GGA calculations, and they have not changed 
their conclusions. However, since the relative energies 
that are used to determine the stability of perovskite sur- 
faces might be influenced by the exchange and correl- 
ation potential, even though DFT+U fails to give better 
results than GGA calculations to predict the phase sta- 
bility of hematite surfaces [19], we still intend to investi- 
gate the effect of DFT+U in later work. The original unit 
cell used to construct the LFO perovskite surface was a 
GdFe0 3 -type orthorhombic unit cell (adapted from Fig- 
ure one in [13]), in which the local magnetic moments 
of Fe are aligned in G-type anti-ferromagnetic order. 
The relaxed lattice constants for a, b, and c in bulk LFO 
correspond to 0.575, 0.559, and 0.792 nm, respectively, 
which are in reasonable agreement with the experimen- 
tal values [20] of 0.558, 0.556, and 0.785 nm. The cutoff 
energies for the wave function and augmentation charge 
density are 25 Ry for the former and 225 Ry for the 
latter. 

We modeled the y/2 x y/2 LFO (001) surface by using 
a repeated slab model. Hamada et al. [10] had already 
shown and we confirmed [13] that one V Q formed in the 
LFO (001) surface promoted the tendency of Pd to 
segregate in bulk. Moreover, we further demonstrated 
that Pd has the strongest tendency to segregate at Fe0 2 - 
terminated surfaces containing V Q s, in comparison with 
three other surfaces, i.e., LaO-terminated surfaces with 
and without V Q s and the perfect Fe0 2 -terminated sur- 
faces. Additionally, Lee et al. [21] calculated a surface 
phase diagram of the LFO (010) surface and argued that 
the LaO-terminated surface could be predicted to be 
stable at lower temperature (773 K), which was in 
agreement with the previous experimental results mea- 
sured by X-ray photoelectron spectra [22,23]. In con- 
trast, the Fe0 2 -terminated surface became dominant at 
high temperatures (> 1,500 K). Therefore, thermal treat- 
ment at high temperature is essential to make Fe0 2 - 
terminated surfaces more stable. We thus examined 
Fe0 2 -terminated surfaces in this work. The atomic 
configuration for a pristine Fe0 2 -terminated surface is 
in Figure 1, which was obtained with visual molecular 
dynamics [24]. Our repeated slab model consisted of 
nine atomic layers, i.e., five Fe0 2 layers and four LaO 
layers. Further, one unit cell contained eight La atoms, 
10 Fe atoms, and 28 O atoms in total. Brillouin-zone 
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Figure 1 Side views of Fe0 2 -terminated surfaces. A vacuum 
region with a thickness of 1 1 A is placed above the top surface. The 
green, brown, and red spheres correspond to La, Fe, and 0, 
respectively. The computational unit cell used in the present 
research is indicated by the black solid line. The bottom two layers 
were fixed during optimization. 



integration was carried out within a Monkhorst-Pack [25] 
scheme using a uniform (4 x 4 x 1) mesh. We inserted a 
vacuum region of 11 A to minimize the interaction 
between two adjacent slabs. We fixed the two bottom layers 
to the bulk coordinates during the geometry optimizations 
and allowed atomic relaxation for the rest of the layers. 

Results and discussion 

Simplified 2D tables that represent the complicated 
atomic configurations of perovskite surfaces have been 
provided in Figure 2 to clarify the discussion. Configura- 
tions with negative formation energies are more stable 
than the reference configuration. One Pd segregating 
from the third Fe0 2 layer to the surface just releases an 
energy of about 0.08 eV [13] (Figure 2 group I (a) and 
(b)) as we demonstrated without V Q s. However, when 
one Pd has already segregated on the topmost site of a 
perfect LFO surface, the additional Pd prefers to stay in- 
side the bulk rather than segregate onto the surface as 
shown in Figures 2 group I (c) to (e). One first has to de- 
termine the positions of V Q s and Pd atoms in studying 
the effect induced by V Q s on the stability of Pd atoms. 
We have to calculate all the possible configurations 
containing V Q s and Pd. Hamada et al. [10] pointed out 
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Figure 2 Simplified 2D tables that represent complicated structures of perovskite surfaces containing Pd„ (n=1 and 2). Groups I to III are 
for the geometries with no V Q , one V Q , and two V Q s, respectively. The atomic configurations for each group, which are schematically represented 
by the table of panel (a), are indicated by the ball and stick model. The uncapping unit cell is indicated by the black line as seen in Figure 1. The 
rows containing Fe (Pd) in each table represent Fe0 2 (Pd0 2 ) layers, and the vertical lines represent O atoms in Fe0 2 (Pd0 2 ) layers. The horizontal lines 
represent O atoms in LaO layers (La atoms are not explicitly shown). The absence of vertical (horizontal) lines means V Q forming at the surface 
(subsurface) site. The calculated difference in energy (in eV) for each panel relative to the total energy of the surface in panel (a) is also listed. 
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that the most stable site for V Q s is the topmost surface 
for pristine LFO and the subsurface (LaO layer) O site 
for Pd located in the first layer of the LFO surface. We 
considered V Q s formed at those two possible sites along 
with various configurations of Pd atoms at the Fe0 2 - 
terminated surface. We set the first configuration in 
panel (a) in group II to the reference state in which one 
Pd atom was located in the first Fe0 2 layer, the second 
Pd atom was in the third Fe0 2 layer, and a V Q was lo- 
cated in the first LaO layer just under the first Pd. The 
positions of the first Pd atom and V Q were found to 
have the most stable configuration. Positive formation 
energies for panels (i) to (m) in group II indicate that 
V Q s that formed on the topmost surface is unstable. 
However, the most stable state was found with a forma- 
tion energy of about -0.57 eV when a V Q was located at 
the subsurface nearly at the center of two Pd atoms, as 
seen in Figure 2 group II (b). However, one of the Pd 



atoms tended to be buried in the second Fe0 2 layer 
(panel (b)) rather than exposed to the vacuum (panel (c) 
in group II), and the energy discrepancy between panels 
(b) and (c) was as large as 0.58 eV. We analyzed the 
projected density of state (PDOS) of the two Pd atoms 
in the V Q -containing surfaces to understand the origin 
for the difference in stability between panels group II (b) 
and (c). All the results are presented in Figure 3. We 
denoted the Pd located at the top-left site in the unit cell 
in Figure 2 group II (a) to (c) as Pd-1 and the other one 
as Pd-2. Where Pd-2 stayed inside the bulk (Figure 2 
group II (a)), the PDOS of Pd-1 looked similar to that in 
Figure five (e) in [13], i.e., a single Pd at the first Fe0 2 
layer with one V Q beneath it. The V 0 beneath Pd-1 re- 
duces hybridization between the Pd state and O p 
state, leading to significant stabilization of the state. 
The degenerated e g states of 4d-orbitals for Pd-2 are sin- 
gly occupied (Figure 3a 2 ). When Pd-2 replaces the Fe 




-4 -3 -2 -1 

E-E F leV 

Figure 3 Calculated projected density of states (PDOS) of two Pd atoms. Panels (aO to (q) are the PDOSs for Pd-1 located at the top-left 
site of Figure 2 group II (a) to (c). Panels (a 2 ) to (c 2 ) represent the PDOSs of Pd-2, which is located at the third Fe0 2 layer (a 2 ), at the subsurface 
(b 2 ), or the first Fe0 2 layer (c 2 ). Positive (negative) values refer to spin-up (spin-down) states. The line through the zero point on the horizontal 
axis represents the Fermi level, (d) is the top view of the calculated isosurface of the partial charge density of the Pd0 2 -terminated surface 
containing V Q in the spin-up state in an energy window of (-0.1, 0.1) eV. The yellow-curved isosurfaces stand for the charge density of 0.6 a.u." 3 . 
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atom at the second Fe0 2 layer (Figure 2b), the change in 
PDOS of Pd-1 seems rather small and similar to that in 
panel (a!). However, the degeneracy of the e g state is 
lifted for Pd-2 because of the missing apical oxygen 
atom, leading to a downward shift in d^ beneath the 
Fermi level, except for a small antibonding state near the 
Fermi level associated with hybridization between the Pd 
d%z and p state of oxygen atom beneath it. The t 2g states 
are also fully occupied in the form of a stable closed 
shell. The degeneracy of the e g state is lifted due to the 
lowering of symmetry at the surface for Pd-2 located at 
the first Fe0 2 layer (Figure 2 group II (c)). However, as 
there is another O at the subsurface, a much stronger 
antibonding Pd state appears near the Fermi level in 
contrast to that in panel (b 2 ). Additionally, the d xy 
state remarkably increases in energy due to increased 
hybridization between the Vc\-d xy and 0-p y/x states, and 
an especially sharp peak emerges at the Fermi level in 
the spin-up state. The Pd d xy state also appears near the 
Fermi level for Pd-1 as shown in panel (ci). The corre- 
sponding partial charge density for the peak at the Fermi 
level has been drawn on the (001) plane in panel (d). 
The spin-up partial charge density exhibits strong anti- 
bonding states in the form of pdjf bonds between Pd 
and O in the energy window from -0.1 to +0.1 eV rela- 
tive to the Fermi energy. As a result, the additional Pd at 
the neighboring surface site is less stable than that at the 
second Fe0 2 layer. 

Since V Q s are more likely to form at the subsurface 
(LaO layer) than the surface in the Pd-containing Fe0 2 - 
terminated surface, we placed another V Q in the same 
LaO layer (Figure 2 group III (a) to (c)). If two V Q s are 
both located at the subsurface, the second Pd atom 
tends to substitute the Fe atom either at the second 
Fe0 2 layer forming a pair of Pd atoms (Figure 2 group 
III (b)) or on the surface forming the Pd0 2 layer 
(Figure 2 group III (c)). The difference in energy 
between these two configurations is less than 0.05 eV. 
Thus, the additional V Q stabilizes the Pd0 2 -layer 
exposed to the vacuum. 

Thus far, we have assumed the existence of V Q - How- 
ever, the concentration of V Q s depends on their forma- 
tion energy. Therefore, we have to verify the stability of 
surfaces containing V Q (s) with different concentrations 
of Pd by taking the formation energy of V Q s into ac- 
count to further strengthen our conclusion. We calcu- 
lated the relative energies of surfaces containing Pd m V Qw 
(m =1 and 2 and n =0, 1, and 2) relative to the perfect 
slab (without V Q s) with Pd inside the bulk of LFO (see 
Figure 2 group I (a)). Note that the systems we have 
discussed here are surfaces with Pd atom(s) and V Q (s) 
located on/near the surface. The relative energies (AE rel ) 
as a function of the chemical potential of oxygen (A^ 0 ) 
are shown in Figure 4. The corresponding geometries 



for the Pd m V Qw -containing surfaces are all included in 
Figure 2. Since two Pd atoms fail to segregate near/at 
the surface adjacently without V Q s, the results for the 
Pd 2 -containing perfect surface excluded from Figure 4. 
The AE rel for each colored line is calculated as: 
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where the first items in Equations 1 to 5 on the right- 
hand side are the total energies of the Pd m V 0w - 
containing (m =1 and 2 and n =0, 1, and 2) surfaces, 
with their subscripts describing their compositions. Efjf 



1500 K 



1000 K 




A^/eV 

Figure 4 Calculated relative energies of five LFO surfaces 
containing Pd m V 0 „. This is with respect to the dissolution phase of 
the LaFe^Pdjft slab as a function of A^/ 0 and oxygen partial 
pressure at high temperatures. 
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represents the total energy of the reference surface that 
contains one solid-solution state of Pd inside the bulk. 
The £ , (F t eo 2 ) 2 /(LaFe03) 8 * n Equations 4 and 5 is the total en- 
ergy of the pristine Fe0 2 -terminated surface. The fi Q is 
the chemical potential of oxygen. The chemical poten- 
tials of oxygen in the LFO bulk and gas phase are equal 
under equilibrium conditions. The fi 0 based on an ideal 
gas approximation is directly connected with the partial 
pressure (p (0 2 )) and temperature (T) by 

<(t» = <(r,/) + *rm(^) , (6) 

= A^(T,p°)+ 1 -kTln(£), (7) 

in which k is the Boltzmann constant and p° is taken to 
be the standard pressure, j/q is the total energy of an 
isolated 0 2 molecule. The A[/q '(T \p 0 ) item is deter- 
mined by using thermodynamic data from the thermo- 
chemical tables [26]. Hence, we can obtain the 
formation energy of V Q (s) based on Equations 2 to 5 by 
subtracting the item in Equation 1. As we assumed that 
the chemical potential of La would be very low in the 
present research, no metallic La or La oxides could pre- 
cipitate on the surface although the La-O bonds were 
broken due to the formation of V Q s in the first LaO 
layer (Figures 2 groups II and III). According to the 
equations, the positive AE rel means the reference surface 
is more stable. 

We can find from Figure 4 that when A^ Q is greater 
than -1.17 eV (point A), no V Q s form on the surface. 
The Pd-segregated surface (Figure 2 group I (b)) is 
slightly more stable than the surface with Pd inside the 
bulk of the perovskite (Figure 2 group I (a)). This indi- 
cates that Pd preferentially stays at the first layer of the 
LFO surface than the bulk position to some extent. One 
V 0 in the surface appears at the subsurface (LaO layer) 
when A/u 0 is lower than -1.17 eV. The surface 
containing Pd 2 V Q is predicted to be stable between 
points A and B, indicating conditions with standard 
pressure at temperatures between 1,000 and 1,500 K. 
Two Pd atoms attract each other in such a surface by 
sharing one V Q in the first LaO layer (Figure 2 group II 
(b)). The PdiVoi-containing surface (Figure 2 group II 
(n)) becomes dominant at A^ Q below -1.67 eV (point B) 
under standard pressure at temperatures over 1,500 K. 
Two V Q s -containing surfaces are predicted to be dra- 
matically unstable compared with the other three sur- 
faces due to the greater formation energy of two V Q s 
under the conditions given in Figure 4. The Pd^c^- 



containing surface (Figure 2 group III (d)) will appear 
under standard pressure at temperatures far above 1,500 
K (the pink line: the critical point is beyond the scale of 
Figure 4). The surface containing Pd 2 V Q2 (Figure 2 
group III (b)) for the blue line is predicted to be unstable 
under any conditions as presented in Figure 4. From 
what we have mentioned above, one V Q can be produced 
at the first LaO layer of the Fe0 2 -terminated surfaces 
with segregated Pd m (m =1 and 2) under reasonable 
working conditions, and such surfaces are predicted to 
be dominantly stable over a wide range of A/u 0 . 

Conclusions 

We investigated what effect oxygen vacancies had on the 
tendency of additional Pd atoms to segregate at the 
LaFe^Pd^Os.y surface, as well as compared the relative 
stability of Fe0 2 -terminated surfaces that contained 
Pd m V Qw versus the oxygen chemical potential, by using 
first-principles theoretical calculations. We pointed out 
that Pd atoms repulse one another without V Q s. How- 
ever, if there are V Q s at the subsurface layer, Pd atoms 
become attractive, forming a pair of Pd atoms while 
sharing one V<> Furthermore, we clarified that the 
Fe0 2 -terminated surface containing Pd m V Q could be 
predicted to become stable over a wide range of oxygen 
chemical potentials below -1.17 eV. Therefore, the 
present results provide support that the calcination of 
precious metals containing catalysts at 1,073 K or high 
temperatures in air during the catalyst preparation step 
leads to the formation of oxygen vacancies near the sur- 
face and then enhances the formation of a LaPdOs^ 
layer in the vicinity of the LaFe0 3 oxide surface. 
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